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ABSTRACT 

We divide some of the observed diffuse interstellar bands (DIBs) into families which 
appear to have spectral structures of single species. Three different methods are applied 
to separate such families, exploring the best approach for future investigations of this 
type. Starting with statistical treatment of the data, we found that statistical methods 
by themselves give insufficient results. Two other methods of data analysis ("Averaging 
EWs" and "Investigating the figures with arranged spectrograms") were found to be 
more useful as tools for finding the spectroscopic families of DIBs. On the basis of 
these methods, we suggest some candidates as "relatives" of 5780 and 5797 bands. 
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1 INTRODUCTION 

Diffuse interstellar bands (see e.g. Herbig 1995), absorption 
structures of interstellar origin still await explanation. The 
identification of the carrier(s) of DIBs is one of the most 
difficult challenges for spectroscopists. 

To the present day, huge amounts of observational data 
on DIBs have been stored by astronomers and published in 
hundreds of papers. Unfortunately, astronomical data still 
do not meet sufficient understanding from the point of view 
of laboratory spectroscopists. 

One tries, in general, to solve the mystery of the car- 
riers) of DIBs on the field of interdisciplinary spectro- 
scopic collaboration between molecular physicists, molecular 
chemists and astronomers. One expects that some progress 
on this field will be possible when all known DIBs (about 
three hundred of them have been discovered in visible light 
region) are divided into families in such a way that only one 
carrier is responsible for all bands belonging to a given fam- 
ily. Such families of bands we call 'spectroscopic' ones, to 
distinguish them from 'characterological' families isolated 
by the other authors (Chlewicki et.al. 1986, Krelowski & 
Walker 1987, Josafatson & Snow 1987). All bands belong- 
ing to the spectroscopic family are, by definition, caused by 
the same carrier. Bands belonging to the characterological 
family have some common characteristics (e.g. all are rather 
narrow), but they may be produced by different carriers. 

To isolate spectroscopic families of bands, first of all 
one has to by-pass in some way the problem of the "noisy 
correlation". In a given spectrogram we have to deal with 
a complicated mixture of interstellar absorption structures. 
This is because the medium between the target star and 



the observer contains various species. For different direc- 
tions (various target stars) we have to deal with different 
column densities of interstellar matter giving contributions 
to the spectra. Intensities of all spectral lines (bands) mea- 
sured in a spectrogram are well correlated with the column 
density of relevant matter and therefore also lines originated 
by different species are mutually correlated. Such correlation 
we call " noisy correlation" . Of course not only differences in 
column densities of relevant interstellar matter may produce 
noisy correlation between DIBs. Other contribution to cor- 
relation of this kind may be given, for example, by mutually 
dependent astrochemical processes. Noisy correlation tells 
us nothing about spectroscopic families of bands. 



The number of recognized DIBs has grown dramatically 
in recent years, primarily due to better quality of observa- 
tional material. More and more weak DIBs (WDIBs) seem to 
appear every time a given spectral region is analysed care- 
fully. One of the authors of this paper (BW) spent a few 
years analysing spectra for new, very weak DIBs. It turned 
out to be of great importance to reexamine the problem of 
DIB families. The aim of this paper is to explore further the 
properties of DIBs in the context of isolating families of the 
structures. We first describe observational material which 
revealed the above mentioned absorption bands. Then we 
present the results of the measuring procedure and describe 
DIB searching methods. In the last section of the paper we 
discuss the problem of some adequate method for separating 
spectroscopic families of DIBs, and we pick out two proba- 
ble "relatives" for 5780 band and four other ones for 5797 
band. 
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Table 1. Measured strengths of DIBs in the McDonald data. HD numbers for target stars are given in the first column. The first row of 
the table contains the names of DIBs (their approximate positions). Individual cells contain subsequently (going from the top): equivalent 
width (EW), roughly estimated measuring error for EW (EW and its error are given in milliangstroms) and line depth, given in promilles 
of the local continuum level. Empty cells in the table indicate to the cases of very bad quality spectrograms. 
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2 OBSERVATIONAL MATERIAL 

All the spectra analysed for the purpose of this paper were 
taken from the archives of Prof. Jacek Krelowski (Astronom- 
ical Center, Nicolaus Copernicus University, Toruh, Poland). 
We used spectra taken with the Canada-France-Hawaii Tele- 
scope (CFHT) (described by Snow and Seab (1991)) cover- 
ing the spectral range 5780-5905 A and spectra acquired at 
the McDonald Observatory with an echelle spectrograph fed 
with the 2.1 m telescope, covering the spectral range from 
5600 to 7000 A (Krelowski and Sneden, 1993). Spectra ac- 
quired with the CFHT have S /N of about 800 and resolution 
R of about 60000. The relatively high quality of the observa- 
tional data used allowed us to measure the equivalent widths 
for many strong and weak DIBs. 



3 MEASUREMENTS 

Firstly, from the McDonald Observatory data, spectrograms 
of the highest possible quality were selected for 25 targets. 
Using the chosen spectra, we measured equivalent widths 
(EW) end depths of 18 prominent DIBs within the spec- 
tral range from 5760 to 6700 A. HD numbers of stars and 
designations of considered DIBs, as well as values of mea- 
sured parameters are given in Table 1. In each cell of the 
table (provided it is not empty) we have (from top to bot- 
tom): EW (in milliangstroms) , error of the measurement of 
the EW, and line depth (in promilles of continuum). In the 
cases of problematic measurements, due to unsatisfactory 
quality of spectra, the cell in the table was left empty. 

We have also selected spectra (McDonald) covering 
spectral range from 5750 to 5860 A for 74 stars and we have 
measured EWs for 11 WDIBs (listed in Table 2) as well 
as for three well known and relatively strong DIBs, namely 
5780, 5797 and 5849. Results of the measurements are given 
in Table 3. Empty cells in the table correspond to situations 
where the measuring procedure was impossible to carry out. 
Errors for EWs in the case of WDIBs are of course large and 
in many cases they correspond to the measurement value. 

Positions of WDIBs given in Table 2 were measured 
in a spectrogram which was averaged over 48 stars (in the 
case of 5760 and 5763 averaging was over 3 stars) observed 
with CFHT. Before performing the averaging procedure, all 
points in individual spectra were shifted along the A-axis a 
value AA (AA/A = const) calculated in such a way that the 
position of the interstellar Nal D2 line became the same as 
in laboratory (A = 5889.95). All measurements were carried 
out by one of the authors (BW) and by using the same mea- 
suring procedure in all cases. The averaged spectrum com- 
pared with the single spectrum of star HD 147165 is shown 
on Fig.l. 



4 SEARCHING THE WDIBS 

The reality of the considered WDIBs and their interstellar 
origin were checked carefully long ago by one of the authors 
(BW) and then it was also claimed by other authors (e.g. 
Krelowski and Sneden 1993, Galazutdinow et.al. 2000). The 
routine tests usually used to prove the interstellar origin of 
the absorption lines were described in the quoted papers in 



detail. Here we focus our attention only on positions and 
rough profile characteristics of WDIBs considered. 

The knowledge of WDIB profiles is of fundamental im- 
portance; firstly, because we need their shapes to measure 
EWs more precisely, and secondly, to make spectroscopic in- 
vestigations easier (mainly the comparisons with laboratory 
spectra of different molecular samples). To study the WDIB 
profiles one needs spectrograms of very high resolutions as 
well as high S/N ratios. Unfortunately, our data are not as 
good as required for such a purpose. To gain higher S/N 
ratio in the spectral range occupied by considered WDIBs 
we have averaged CFHT spectra over different stars. Aver- 
aged spectra obtained are displayed on Figure 1. Using these 
spectra we have measured positions of the WDIBs (Table 
2). Measurements of the EWs for WDIBs were carried out 
keeping in mind their profile characteristics, visible in Figure 
1. Profile characteristics visible in averaged spectrum were 
helpfull during EW measuring procedure for WDIBs in in- 
dividual spectra with unresolved left and right limits of the 
profile. 



5 STATISTICAL ANALYSIS FOR STRONG 
DIBS 

As far as relatively strong bands are concerned we have anal- 
ysed statistically mutual dependence between EWs, and sep- 
arately between depths, for 18 prominent DIBs (data from 
Table 1). Because all considered bands are very far from 
being saturated, the members of one spectroscopic family 
should be linearly correlated. 

One of the most important tests giving information 
about linear dependence between two variables is the test for 
the existence of linear correlation between these variables. 
Therefore, we obtained coefficients of linear correlation r, 
for all considered pairs of bands. 

Results of the study of EWs and depths for 18 promi- 
nent DIBs are given in Table 4. Preliminary analysis show 
that all prominent DIBs are mutually dependent. Thus, dur- 
ing our search for possible spectroscopic families, we should 
concentrate on the cases of extremely high correlations. 
Marks in Table 4 were put only for these cases of correlation 
for which the value of statistic t (see Brandt 1970): 



(where R is a value of the parameter of the linear corre- 
lation obtained during estimation) distributed according to 
Student distribution with n — 2 degrees of freedom, was twice 
as much as that obtained for the confidence level 0.995. 

Theoretically on the basis of results contained in Table 
4 we would be able to select some DIB families. However, 
by inspection of Table 4, we find for example, that the 6196 
band has in its family also the following bands: 6203, 6270, 
6426, 6614, 6699, 5780, 5797. But the 5780 and 5797 bands 
do not belong to the same spectroscopic family as was shown 
by Krelowski and Westerlund (1988). This indicates clearly 
that DIBs could be quite well correlated even if they do not 
belong to the same spectroscopic family, and that we cannot 
rely only on a formal statistical approach when looking for 
such families of DIBs. It is very probable that abundances 
of different agents of DIBs are relatively well correlated only 
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Figure 1. Averaged (over different stars) CFHT spectrograms. Vertical arrows indicate positions of studied bands. The shorter piece of 
spectrogram is a result of averaging over 3 stars (but these spectrograms were of the highest signal to noisy (S/N) ratios) whereas the 
longer one gives an average over 48 spectrograms of moderate S/N ratios. Spectrograms used in the averaging procedure look like this 
one for HD147165, given here for comparison. 
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Table 2. The list of names and exact positions for WDIBs. Positions were measured in averaged spectra (Fig.l) and their errors, in 
most cases, should be less than 0.1 A. 



due to properties of interstellar matter. Physical conditions 
which allow one molecule to be formed within interstellar 
clouds simultaneously may give a chance for the other ones 
to also originate. In such cases certain noisy correlation has 
to appear. 



6 LOOKING FOR CORRELATIONS 

BETWEEN WDIBS AND 5780 AND 5797 

The extremely low EWs of WDIBs (in most cases less than 
10 mA) make their measurements very difficult. Avoiding 
large errors in the measured parameters is impossible. In 
such a case one looks for an adequate method to find any 
possible correlations between individual WDIBs and cho- 
sen strong DIBs. As it was shown in section 5, statistical 
correlations do not always give us entirely satisfying results 
even if the considered bands are relatively strong. As al- 
ready mentioned, 5780 and 5797 bands do not belong to 
the same spectroscopic family (Krelowski and Westerlund 
1988). Each of these strong bands may therefore play a role 
of a representative band for its own family and may have its 
"relatives" among WDIBs. 

6.1 Statistical analysis of the intercorrelations 
between WDIBs and prominent DIBs 

We have analysed statistically mutual dependence between 
strong and weak DIBs (data from Table 3). In some cases 
we obtained from measurements value zero or had a lack of 
measurements for WDIBs. We excluded all such cases from 
the analysis. Computed coefficients of linear correlation be- 
tween EWs of each of the three strong bands (5780, 5797, 
5850) and EWs of WDIBs are presented in Table 5. In all 
cases (with only one exception 5780 versus 5769) the param- 
eter of linear correlation r is greater than zero at a high sig- 
nificance level (for level of confidence 0.95 the critical value 
of r = 0.23). This means that EWs of each of the strong 
bands are mutually dependent (in the statistical sense) and 
correlated. However, it should be noted that in some cases, 
even if the value of the parameter r is high enough accord- 
ing to our criteria, the \ 2 test suggests that a more complex 
dependence, rather than the linear one, may be relevant. 
From the detailed analysis we obtained: (i) WDIBs 5760, 
5763, 5809, 5819 are linearly-dependent with all considered 
strong bands, as was checked using the \ 2 test, (ii) in the 
case of strong bands 5797 and 5850 we can also accept the 



hypothesis of the linear regression fit for the 5793 and 5829 
WDIBs, contrary to the strong band 5780 for which such hy- 
pothesis does not work, (iii) for the other bands (i.e. 5766, 
5769, 5773, 5776 and 5795) the hypothesis of linear depen- 
dence with all considered strong bands should be rejected. 

The result in point (i) clearly shows that some of 
the WDIBs are linearly well correlated with different DIBs 
(5780, 5797) which definitively do not belong to the same 
spectroscopic family. Therefore, the statistical approach to 
our data is not sufficiently strong by itself to discriminate the 
membership of WDIBs into families. More detailed analysis, 
taking into account other methods, is therefore required. 



6.2 Averaging EWs method 

To check mutual dependence between WDIBs and 5780 and 
5797 bands, we have started from the following procedure. 
From Table 3 we extracted a sample with EW(5797) re- 
stricted to a very narrow range of values. From such a sam- 
ple we separated two subsamples: the first one with rela- 
tively low, and the second one with relatively high values 
of EW(5780). In these subsamples we calculated mean val- 
ues of EW(5780), EW(5797) and EWs for WDIBs. Both 
subsamples, as well as results of averaging, are displayed in 
Table 6. WDIBs correlated with 5780 should follow changes 
in EW(5780), contrary to those which are not correlated. 

In a quite similar way we have extracted a sample of 
stars with almost constant EW(5780) and separated sub- 
samples with opposite values of EW(5797) - Table 7. Such 
a simple procedure allowed us to divide WDIBs into three 
groups; (i) well correlated with 5780 (5776, 5795), (ii) well 
correlated with 5797 (5793, 5819, 5829) and (iii) correlated 
with neither 5780 nor 5797 (5769, 5773, 5809). Three other 
cases (5760, 5763, 5766) are rather complicated because they 
seem to be correlated with 5797 as well as with 5780. The 
main selection criterion here was a gradient of the average 
value for WDIBs. If this gradient, for a considered WDIB, 
is evidently similar to the gradient which is typical for the 
leading band (in our case this gradient is very close to 0.5) 
we included this WDIB into the family. 

Please note that the term "correlation" used above does 
not mean the same as correlation in the statistical sense. 
Statistically, all considered bands are quite well mutually 
correlated. 
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Table 3. Measured EWs (in mA) for WDIBs listed in Table 2 and for three prominent bands - 5780, 5797 and 5850. Empty cells 
correspond to cases of quality of spectrograms too bad to find any weak band. Zero values were put in for non-remarkable bands in 
spectrograms of quite good quality. 
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Table 4. Line-line correlations between prominent DIBs listed in Table 1. The sign "+" means that EWs of considered lines are well 
correlated, with the value of coefficient of linear correlation, r, greater than 0.77. When this coefficient is greater than 0.9 the sign "*" 
is used. The sign "X" stands for good correlation (r > 0.79) between the depths of bands. 
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Table 5. The correlation table between strong DIBs (listed in the first row of the table) and WDIBs (listed in the first column of the 
table). Formal coefficients of linear correlation were calculated on the basis of the data taken from Table 3. 



6.3 Investigating the figures with arranged 
spectrograms 



Another way to divide WDIBs into the three above men- 
tioned groups consists in visual inspection of spectrograms. 
If we compose sequences of spectrograms in such a way that 
EW(5780), for example, is constant and EW(5797) grad- 
ually increases (Fig. 2) we are able, after visual inspection 
of the picture, to indicate WDIBs which follow, and those 
ones which do not follow, the behaviour of 5797. Figure 3 



shows the case of nearly constant EW(5797) and changing 
EW(5780). 

This inspection procedure of Figures 2 and 3 confirms 
results achieved with the earlier described averaging proce- 
dure. Furthermore, one can, when studying Figures 2 and 3, 
make the following conclusions: 

- the 5850 band, which corresponds to relatively strong 
phenomena, is very well correlated with the 5797 band and 
not (at all) with the 5780 one (this was claimed already by 
Krelowski et al, 1993), 
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Table 6. Two subsamplcs of stars chosen in such a way that average EW(5797) is almost constant in both subsamples, contrary to 
EW(5780) which in the first subsample is almost twice as large as in the second one. The names of bands are listed in the first row, 
whereas names of target stars used are given in the first column of the table. In both subsamples for each band the mean value of 
measured EWs was counted. The sign "-" stands for the lack of the data and "+" indicates interesting examples of correlation (see text). 
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Table 7. The same as in Table 6 but here the mean E W(5780) is rather stable whereas E W(5797) drops almost in half when going from 
the first subsample to the second one. 
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Figure 2. CFHT (the longest pieces) and McDonald's (shorter, individual orders of echelle spectrograph) spectrograms arranged in 
a way that strength of 5797 band gradually increases whereas intensity of 5780 band remains constant. 
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- the 5763 band follows the behaviour of the band 5797 
and does not follow (at all) the 5780 band, 

- although Figures 2 and 3 are not as good as required 
for studying the 5760 band it seems that this WDIB is cor- 
related rather with 5797 than with the 5780 band, 

- the 5769 band also seems to be well correlated with 
the 5797 band. 

It is worthly of note that the observations we are dealing 
with here were not planned for searching for spectroscopic 
families. One may, however, perform much more adequate 
observations, directed to the subject. In the case of better 
data, similar analysis would give much more reliable results 
and lead directly to separation of true spectroscopic families 
of bands. 



7 DISCUSSION 

As mentioned in the introductory section, the main obsta- 
cle to isolating spectroscopic families of DIBs is the noisy 
correlation. Due to a very high level (as the analysed data 
reveals) of noisy correlation in the considered data, the abil- 
ity of statistical method to isolate spectroscopic families are 
very limited. 

Tight linear correlation, expected between members of 
the same spectroscopic family, is efectively hidden by noisy 
correlation and by measuring errors. This is evident when we 
study results of statistical analysis described in the previous 
section. Looking at Table 5 only, we are not able to indi- 
cate which WDIB belongs to the same spectroscopic family 
as, for example, the band 5797. On the other hand, when 
considering Table 4, we would be inclined to isolate, e.g, the 
family of strong lines: 6196, 6203, 6270, 6426, 6614, 6699, 
5780, 5797; and this would be a mistake, since 5780 and 
5797 belong to different families, as mentioned in section 5. 

Taking into account that statistical analysis requires 
plenty of usable data and gives insufficient (for solving our 
problem by itself) results, it is therefore not recommended 
as an appropriate tool for isolating spectroscopic families. 
[However, the statistical approach may be useful to distin- 
guish linearly correlated bands from these ones which are 
correlated non-linearly. Non-linearly correlated DIBs should 
belong to different spectroscopic families. Also multidimen- 
sional statistical analysis could be useful in this case. Further 
study of this problem will be the subject of a separate paper 
(Godlowski & Wszolek, in preparation)]. 

Methods described in subsections 6.2 and 6.3 are much 
more appropriate than the formal statistical aproach. They 
need a fewer number of spectrograms and much less time 
for making EW (or depth) measurements. These methods 
seem also to give valuable results. Using these methods we 
performed preliminary separation of two presumable spec- 
troscopic families: (i) 5780, 5776 and 5795, and (ii) 5793, 
5797, 5819, 5829 and 5850. 

Most probably, designated families are not complete 
yet. One cannot exclude also the possibility that we made 
wrong indications. In the case of almost constant ratios be- 
tween column densities of various DIB carriers in interstellar 
clouds, we have a chance to get results quite similar those 
of the case when we have to deal with few spectral lines of 
the same carrier. Further investigation, based on better data 



samples and involving other spectral ranges, is necessary to 
isolate true spectroscopic families of bands. 
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